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Abstract 

We study a triple singular limit for the scaled barotropic Navier-Stokes system modeling the 
motion of a rotating, compressible, and viscous fluid, where the Mach and Rossby numbers are 
proportional to a small parameter e, while the Reynolds number becomes infinite for e — t- 0. If 
the fluid is confined to an infinite slab bounded above and below by two parallel planes, the 
limit behavior is identified as a purely horizontal motion of an incompressible inviscid fluid, the 
evolution of which is described by an analogue of the Euler system. 
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1 Introduction 

There have been several types of singular limits of the Navier-Stokes system studied in the recent 
literature, among them a few devoted to the effect of rotation, see Chemin et al. [1]. In this paper, 
we consider the scaled compressible Navier-Stokes system describing the time evolution of the density 
Q = g(t, x) and the velocity u = u(t, x) of a compressible viscous and rotating fluid: 

dtg + div-,{gu) = 0, (1.1) 

^((^u) + div^(^u (g) u) + r^g{oj X u) + z^VxP{g) = r^div^S + gVxG, (1.2) 

Ro Ma Re 

where § is the viscous stress, here given by Newton's rheological law, 

u + V^u - ^div^ul") + r/div^.u| , (1.3) 

P = p{q) is the pressure, fi > and rj > are the viscosity coefficients, u = [0, 0, 1] is the axis of 
rotation, and VxG represents a conservative force imposed on the system, say, by the gravitational 
potential G of an object placed outside the fluid domain, see the survey of Klein |21] . 
The scaled system contains several characteristic numbers: 

• Ro Rossby number 

• Ma Mach number 
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Re 



Reynolds number 



The following are examples of singular limits considered in numerous studies: 

• The low Mach number limit: The Mach number is the ratio of the characteristic speed of 
the fluid divided on the speed of sound. In the low Mach number limit, the fluid flow becomes 
incompressible, the density distribution is constant and the velocity field solenoidal, see Ebin 
[7], Klainerman and Majda [12], Lions and Masmoudi [22], Masmoudi [2H], among others. 

• The low Rossby number limit: Low Rossby number corresponds to fast rotation. As 
observed by many authors, the highly rotating fluids become planar (two-dimensional). Ac- 
cordingly, the fast rotation has a regularizing effect, see Babin et al. [I], [2], Chemin et al. 

• The high Reynolds number limit: In the high Reynolds number limit, the viscosity of 
the fluid becomes negligible. Consequently, solutions of the Navier-Stokes system tend to the 
solutions of the Euler system, see Clopeau, Mikelic, Robert [S], Masmoudi [23], [2S], [2S], 
Swann [21], among others. The inviscid limits include the difficulties related to the boundary 
behavior of the fluid and a proper choice of boundary conditions, see Kato [17] , Kelliher [18] , 
Sammartino and Caflisch [29], [30], Temam and Wang [32], [33] . 

The effects described above may act simultaneously. The incompressible inviscid limit was inves- 
tigated by Masmoudi [25], for viscous rotating fluids see Masmoudi [21], Ngo [2Z], the compressible 
rotating fluids were discussed in [9], [10]. In this paper, we address the problem of the triple limit 
for Ma = Ro = e — )■ 0, while Re = Re(e) — oo as e — > 0. In agreement with the previous discussion, 
the fluid flow is expected to become (i) incompressible, (ii) planar (2D), and (iii) inviscid and as 
such described by a variant of the 2D incompressible Euler system that is known to possess global-in- 
time solutions for any regular initial data. Note that the action of volume forces in the momentum 
equation (11. 2p is represented solely by the potential G, notably the effect of the centrifugal force 
is neglected. This is a standard simplification adopted, for instance, in models of atmosphere or 
astrophysics, see Jones et al. [TJ], [IS] , Klein [20]- On the other hand, although the centrifugal force 
is counterbalanced by the gravity in many real- world applications (see Durran [6]), it is proportional 
to 1/e^ under the present scaling, and its far- field impact may change the limit problem dramatically, 
see [9]. 

Similarly to [H], [ID], we consider the problem (11.11 - II. 3p in an infinite slab Q = x (0,1), 
denoting the horizontal components of a vector field v by v/j, v = [v/i,?;3], where the velocity field u 
satisfies the complete slip boundary conditions 

u - n = uslen = 0, [§(V,u)n]tan = [^2,3, -5i,3, 0] |an = 0, (1.4) 
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where n = [0, 0, ±1] is the outer normal vector. Such a choice of boundary behavior prevents the 
flow from creating a viscous boundary layer - the up to now unsurmountable difficulty of the inviscid 
limits, see Kato [17], Temam and Wang [33]. As a matter of fact, replacing (11. 4p by the more standard 
no-slip boundary condition would drive the fluid to the trivial state u = in the asymptotic limit 
unless we impose anisotropic viscosity, see Bresch, Desjardins and Gerard- Varet [3], Chemin et al. 

Our approach is based on the relative entropy inequality (cf. [TT], Masmoudi [2S], Wang and 
Jiang [35) applied in the framework of weak solutions to the Navier-Stokes system (11.11 - [L3|) . We 
consider the ill-prepared initial data: 

f?(0,-) = Qo,e = Q + eg^l f?g ^ in L'{n), (1.5) 
u(0,-) = uo,e^uoinL2(fi;/23)^ (1.6) 

where > is the anticipated constant limit density enforced by the incompressible limit. Accord- 
ingly, the resulting problem is supplemented by the far fleld conditions 

Q ^'Q, u — 7- as \xh\ — )■ oo. (1.7) 

Supposing we already know that, in some sense, 

(I) Qe-Q ^ 
Qe = > 1^ ^ V 

we may (formally) check that q = q{xh), v = [v/i(a;/i), 0] are interrelated through the diagnostic 
equation 

X V + ^^V^g = 0, (1.8) 
Q 

and q satisfles 

i^'" - ^)') ^ ■ (^^ - vm') = °- 

Here and hereafter, the subscript h indicates the restriction of the standard differential operators to 
the horizontal variables, for instance, V^,/ = [d^^f ,8x2!], div/^v = d^^hi + dx2h2, = div/jV/i, etc. 
Note that _ 

P\q) 

whence r can be viewed as a kind of stream function, while the system (II. 8p . (I1.9p possesses the 
same structure as the 2D Euler equations. In particular, we expect the solutions of (II. 8p . (II. 9p to 
be as regular as the initial data and to exist globally in time. Equation (II. 9p arises in the theory of 
quasi-geostrophic flows, see Zeitlin [35l Chapters 1,2]. 
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One of the major stumbling blocks in the analysis of the singular limit is the presence of rapidly 
oscillating Rossby-acoustic waves. Their behavior is described by means of a hyperbolic system 



edts + ^diVj^V = 0, 



(1.10) 



edtY + w X V + 



e 



) 



(1.11) 



where 



Qe- Q 



V ^ V, = u,. 



e 



The bulk of the paper is devoted to the dispersive estimates for the problem fll.lOp . f ll.lip . In 
particular, we use the recent results of Guo, Peng, and Wang [H] on the asymptotic behavior of the 
abstract group of operators 



where $ is a function with specific properties. In particular, we establish — L°° decay estimates 
for the solutions of (11.101) . (II. lip in the frequency domain bounded away from zero. 

The paper is organized as follows. In Section [2], we introduce the standard definition of finite- 
energy weak solutions to the scaled system (11.11 - 11.60 and formulate the main result. Section [3] 
contains the relative entropy inequality, together with the uniform bounds on the family of solutions 
of the scaled system. Section H] represents the bulk of the paper. Using the abstract result of Guo et 
al. [H] we establish the L^ — L°° estimates for the acoustic- Rossby waves. Such a result for the system 
(ll.lOp . (II. lip may be of independent interest and represents an analogue of the standard Strichartz 
estimates for the wave and Schroedinger equations. In particular, we extend the smoothing estimates 
established in [TU] and obtain the necessary tool to attack the inviscid limit in Section 13 

2 Preliminaries, main results 

In order to fix ideas and to simplify presentation, we suppose, without loss of generality, that ^ = 1. 
In addition, we assume that the pressure p G C[0, oo) fl C^(0, oo) satisfies 




p(0) = 0, p'{q) > for all g > 0, lim 




Poo > for a certain 7 > -. 



(2.1) 



Morever, again for the sake of simplicity, we suppose 



p'Cq) = 1- 



(2.2) 
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Finally, given our choice of the complete slip boundary conditions f ll.4p . it is convenient to replace 
the set n = R'^ X [0, 1] by 

meaning we suppose that all quantities are 2-periodic with respect to the vertical variable 0:3. More- 
over, in accordance with fll.4p . we assume that 

0{t, Xh, -X3) = 0{t, Xh, X3), Uh{t, X, -X3) = Uh{t, X, X3), us{t, X, -X3) = -usit, X, -X3), (2.3) 

and 

G{xh, -X3) = G{xh, X3) (2.4) 

for all t G (0,T), Xh € R^, X3 G [— 1, Such a formulation, completely equivalent to (11. 4p . 

was proposed by Ebin [S]. 

Setting Ma = Ro = e, = /i^ \ we say that f), u is a finite energy weak solution to the scaled 
Navier-Stokes system (11.11 - II. 7p if: 

• The density g is a. non-negative function such that 

{g-l)eL^{0,T; [L' + L^){Q))- 

the velocity u belongs to the space L^(0,T; W^''^{Q)). Moreover, in accordance with our con- 
vention, both g and u satisfy the symmetry condition (12.31) . 

• The equation of continuity (II. ip holds in the weak sense: 

/ / {Qdt^p + gu ■ V^ip) dx dt = - [ go,e^{0, ■) dx (2.5) 
JO Jn Jn 

for any ip E C^{[0,T) x fi). 

• Similarly, the momentum equation is replaced by a family of integral identities 

(^g\i ■ dtif + gu^u: Vx-v? - x u) ■ (p + ^p{g)divx(p^ dx dt (2.6) 

= / / Se(Vx.u) : V^ip dx dt- f f gV^G ■ (p dx dt - f ^io,£Uo,e ■ <piO, ■) dx 
Jo Jn Jo Jn Jn 

for any ip E G^{[0,T) x n), with 

§,(Vxu) = fi, (v^u + V^u - ^div^-ul) , fie\0. (2.7) 
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The energy inequality 
■1 , ,o 1 



/ g\u\' + -{H{g)-H'{l){g-l)-H{l)) {r, ■) dx + T f S,{V,u) : V.u dx dt (2.8) 
Jn 12 J Jo Jn 

< f \^g,^^\noJ' + \{H{go.e)-H'{l){go,e-l)-H{l))\ dx + T f gV.G-udxdt 



holds for a. a. r G [0,T], where we have set 



Jl 



(2.9) 



2.1 Limit system 

Under the convention (12.11) . (12.21) . the expected hmit problem reads 

w X V + V^g = 0, V = [vft,(x/i), 0], q = q{xh), 



dt[Ahq-q)+y^h-VhiAhq)=0, 
supplemented with the initial condition 

g(0, ■) = go- 

Note that (12. lip can be written as 

dt^Ahq -q)+ V/, • V/,(A/,g - q) = 0; 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



whence the problem enjoys strong similarity with the standard Euler system. In particular, we may 
use the abstract theory of Oliver [28l Theorem 3] to obtain the following result: 

Proposition 2.1 Suppose that 

qo e W'^'^R^) form> 4. 
Then the problem ^2.12^ . 1^2. 13\} admits a solution q, unique in the class 

q e C{[0,T];W"''^{R')nC\[0,T];W"'-^'\R^)). 



It is worth noting that, similarly to the 2D-Euler system, the solution q can be constructed 
lobally in time. 
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2.2 Main result 

Having collected all the necessary preliminary material we are in a position to state the main result 
of the present paper. 

Theorem 2.1 Let the pressure p satisfy the hypotheses Ii2. Ii2.2\) . Suppose that the initial data 
Qo,e, ^o,e belong to the symmetry class l{2.3\) and are given through U.5\) . U.6\) . where 

{Qoi}e>o bounded in n ^[,^] ^ g'^^ m ^^(fi), (2.14) 

{uo,e}e>o bounded in L^{Q;R^), Uq^^ — Uq in L^(f2;i?^), (2.15) 

where 

G IV-''2(fi), uo G iy'2(fi;i?^), m > 3. (2.16) 

Let 

V^G G L°° n L^(fi; R^) for a certain 1 < r < 2, (2.17) 

and satisfy ^2.4^ . 

Furthermore, let go = Qoi^h) be the unique solution of the elliptic problem 

- Ahqo + qo= t curU[uo],. dxs + /' oi^ dxg m W^^\R^). (2.18) 
JO Jo 

Finally, let [QsiVLs] be a weak solution of the scaled Navier-Stokes system U.l\ - \T7^ in the sense 
specified above. 
Then 

weakly-(*) m L°°(0,T; (L^ + L^)(n)) 



> q 



(strongly) in Li^{{0,T) x il), 
weakly-(*) m L'^{0,T; L^{n; R^)), 



strongly m Ll^{{0,T) x R^) 
where [g, v] is the (unique) solution of the problem ^2.1(A - 12.12\} . 

The rest of the paper is devoted to the proof of Theorem 12.11 
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3 Relative entropy, uniform bounds 



We start by introducing the relative entropy functional for the compressible Navier-Stokes system 
identified in [11], [12], Germain [13]. Set 



J ^ _ ^ s 
where the function was defined in (12.91) . 

3.1 Relative entropy inequality 



r)-H(r)](T,- 



dx, 



(3.1] 



It can be shown that any finite energy weak solution [^£,U£] of the Navier-Stokes system (11.11 - [L7|) 
satisfies the relative entropy inequality: 

£, u, I r, U) (r) + [ (§e(V,.u,) - S,(V,U)) : (V,.u, - V,.U) da; dt < (3.2) 

^.(f?o,.,Uo,. I r(0,-),U(0,-)) 

+ r f Qe {dtV + u, ■ V,U) ■ (U - u,) dx dt 

Jo Jn 

+ r f S,(V,U) : V,(U - u,) dxdt + - r f g,{u x u,) • (U - u,) da; dt 
Jo Jn e Jo Jn 

diY^\j{p{gs) -p(r)) da; dt 



Jo Jn 



dx dt 



Jo Jn 



(r - ge)dtH'{r) + V,H'{r) ■ (rU - 

- r [ g.V.G ■ (U - u,) da; dt 
Jo Jn 

for all (smooth) functions r, U such that 

r > 0, (r - 1) G C~([0,T] x H), U G C^{[0,T] x fi). 



(3.3) 



see [TT]. Clearly, the class of admissible "test" functions (13. 3p can be considerably extended by 
means of a density argument. Note that the relative entropy inequality (13. 2p reduces to the energy 
inequality ( 12. Sp provided we take r = 1, U = 0. 
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3.2 Uniform bounds 

Before deriving the available uniform bounds on the family of solutions [f?e,U£] it seems convenient 
to introduce the essential and residual component of any function h\ 



hess = x{Qe)h, X ^ C*^(0, oo), < X < I5 X = ^ on an open interval contaning ^ = 1, 

/ires = (1 - XiQe))K 

cf. [9], m- 

The uniform bounds are derived from the energy inequality ( 12. 8p (the relative entropy inquality 
(13. 2 p with r = 1, U = 0). Since the initial data satisfy the hypotheses (I2.14p . (I2.15p . the integral on 
the left-hand side of (12. 8 p remains bounded uniformly for e — 0. Accordingly, we get 



ess sup \\y/Q^U^\\^2(^.R3) < c, 
iG{0,T) ^ ' ' 



ess sup 

ie{0,T) 











e . 


ess 



< c, 



L2(n) 



ess sup ||[^£]res||l7(Q) +ess sup \\[l]res\\lj^n) <e^c, 



te{o,T) 



and 



He 



T 



10 Jn 

see [ini Section 2]. We remark that 



t6(0,T) 

, 2 

Va;Ue + Vj^ut div^jU^I 



da; dt < c, 



(3.4) 

(3.5) 
(3.6) 

(3.7) 



g,V,G ■ u, dx| < c (&|u,|2 + {[g,]l^ + 1) \V,G\^) dx, 

where, thanks to the hypothesis (12.170 . the right integral can be "absorbed" by means of a Gronwall- 
type argument. 

Finally, by virtue of Korn's inequality, the relation (13. 7p implies 



T 



Jn 



iVrUpl dx dt < c. 



(3.8) 



Note that all the above bounds depend, in general, on T. 
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3.3 Convergence, part I 

It follows immediately from (13.41 - I3.6p that 

^(1) ^ _^ ^(1) in L~(0,T; (L' + U){n)), r = min{2,7}, (3.9) 

^/^u, ^ u weakly-(*) in L°°(0, T; ^^(fi; R^)) (3.10) 
at least for suitable subsequences. In particular, 

1 in L°°(0,T; (L2 + L")(fi)), r = min{2,7}, (3.11) 

and 

^,u, ^ u weakly-(*) in L^{0,T; {L^ + L^^^^^+^^il; R^)). (3.12) 
Finally, letting e — in the equation of continuity (12. 5 p we deduce 

diVa;U = 0, (3.13) 

and, multiplying the momentum equation by e, we get the diagnostic equation 

w X u + V^^(^) = 0, (3.14) 

where both relations are to be understood in the sense of distributions. It is easy to check that (I3.14p 
imposes the following restrictions: 

g^^^ independent of X3, meaning g^^^ = g^^\xh), (3.15) 

u = u{xh), div^u = divhUh = 0. (3.16) 
Finally, since u belongs to the symmetry class (12.30 . 

M3 = 0, u= [u;„0]. (3.17) 

4 Dispersive estimates 

As already pointed out in the introduction, the heart of the paper are dispersive estimates for the 
acoustic-Rossby waves, the propagation of which is governed by the system 

edts + div^V = 0, s(0, ■) = So, (4.1) 

edtV + X V + V,.s = 0, V(0, ■) = Vq. (4.2) 
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4.1 The wave propagator 

Consider the operator 

B: 



s 




V 


I-)- 



diVj^V 
a; X V + VxS 



defined on the space L^(f2) x L'^{VL\B?). The operator B is skew symmetric, with the domain of 
definition 

V[B] = {[r,V] I r G W^^\n),V G ^^(fi; i?^), div^V G L^(fi)}. 
Next, we introduce the null space N'lB], 

J\f{B) = {[g, v] I q = q{xh), q G W^'\R^), v = [^rh{xh),V3], divh^h = 0, d^.v^ = 0, (4.3) 

w X V + V^g = o|. 

We remark that f 3 = as soon as [g, v] G A/'[;B] belongs to the symmetry class f l2.3p . 
4.1.1 Projection onto N'iB) 

Given a couple of functions [r, U] defined in Q, we want to compute the orthogonal projection 

V : L\n) X L\n-R^) ^Af{B), P[r,U] = [g,v]. 



In addition, we assume the [r, U] belongs to the symmetry class fl2.3p . 

To begin, we project [r, U] onto the space of functions depending only on the horizontal variable 
Xh, meaning we take 

r{xh) = \ I r{xh,X3) dxs = [ r{xh,X3) dxg, IJ h{xh) = [ lJh{xh,Xs) dxg, U3 = 0. 



2 J-i Jo 
Our goal is to minimize the functional 







under the constraint 
We have 



V] ^ \\q - f||i2(R2) + ||V^ - U/,||i2(R2.R2) 

w X V + Vx-g = 0. 



~ll2 



I I 
\1 ~ ^llL2(ij2) + 



+ ^h-Vh 
(w X v) — (w X U) 

12 



L2(/J2.^2) 
2 



2 



= Ik - ^IIl2(r2) + V/ig+(a;xU) 

Thus the associated Euler-Lagrange equation for the minimization problem give rise to 

- AhQ + q= f cmlhUh dx^ + [ r dxa, (4.4) 
JO Jo 

and 

V = [1^1,^2], Vi = -d^^q, V2 = d^,q, (4.5) 
cf. the initial data decomposition fl2.18p . 

4.2 Spectral analysis and dispersive estimates 

We employ the methods of Fourier analysis in order to derive dispersive estimates for solutions of 
the system (14. ip . (14. 2p . Formally, the solutions of (14.11) . (14. 2 p may be written in the form 

[s,Y]{t)=exp(i-[iB])[so,Vo], (4.6) 



e 

where iB is a selfadjoint operator in L^(fi) x L^(f2; R^). 

Accordingly, we deduce that the solution operator generates a group of isometrics in the L^-norm, 
specifically 

Mt, -Winn) + l|V(t, ■)lli2(f,.«3) = ||so||i2(j,) + ||Vo||i2(f,.^3) for any t E R. (4.7) 
Moreover, as the problem is linear, we obtain 

\\s{t, ■)\\wm,2^Q-^ + ||V(t, ■)llw"".2(Q;R3) = ||so||vF™>2(Q) + ||Vo||^,n,2(Q.^3) for auy t e R, m = 0, 1, . . . 

(4.8) 

4.2.1 Fourier representation 

For each function g G L'^{Q), we introduce its Fourier representation 

where ^ ^ 

^(^' k) = ~^ / exp (-i^ ■ Xh) g{xh, xs) dxh exp (-iA;x3) dxs. 
V2 J~i Jb? 

We have 

g{xh, X3) = ^ ^^Xxy, [9{i^ k)] exp (ikx^) , 
where the symbol J^x^^e, denotes the standard Fourier transform on R^. 
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4.2.2 Solutions in the Fourier variables 

The problem fl4.ip . f l4.2p expressed in terms of the Fourier variables reads: 



dt 



whence 









+ i 









6 6 A; 

6 i 

6 -i 

A; 



s{t,^,k) 
V2it,^,k) 



t 



0, 



" siO,^,k) ' 




So{^,k) 






. Vo(e,fc) . 



(4.9) 



[s{t, e, fc), V(t, e, A;)] = exp ( k) ) [so(e, A:), Vo(e, A;)], 



(4.10) 



with the symmetric matrix 



A{C,k) 



^1 ^2 k 

^1 i 

6 -i 

A; 



It is a routine matter to check that the symmetric matrix A{^, k) possesses four eigenvalues: 

1/2 



Ai(e,A:) 

m.k) 



l + l?P+fc2 + ^(l+|g|2+A;2)2_4fc2 



l+|gP+fc^-v/(l+|gP+fc2)^-4fc2 



A2(e,A;) = -Ai(e,A;), 



1/2 



(4.11) 



, A4(e,A;) = -A3(e,A;). 



Note that A3(^,0) = A4(^,0) = 0, which corresponds to the non-trivial kernel of the operator B 
discussed in Section 14.1.11 Consequently, diagonalizing the matrix A., we may rewrite (14.10p in the 
form 

[s(t,e,A;),V(t,e,A;)] (4.12) 
/ \ XiiLk) 1\ 



Q*(e,^)exp 



V 



Ai({,A;) 

-Ai(^,A:) 

A3(e,A;) 

-A3(e,A;) 



Q{iMso{Lk),Yo{L 



for a suitable matrix Q. 
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4.2.3 Decay estimates 

For each fixed k, the solution operators introduced in fl4.12p may be viewed as 

exp (-itAj(e, k)) = exp (-itAj(|e|, k)) ^ exp (^-itAj (^V^^, k) ) , j = 1, . . . , 4. 

In particular, the eigenvalues are smooth functions of |^| ~ ^/—Ah on the open interval (0, oo), and, 
as can be checked by direct computation, Ai(|^|,/c) is strictly increasing in |^| for any fixed k & Z, 
while Xsd^lyk) is strictly decreasing whenever k 0. Consequently, we can use the result of Guo, 
Peng, and Wang [T^ Theorem 1 (a)] to obtain the decay estimate 



exp -itAj A/-Ah, k ^ J -Ah [v 



J = 1,2, keZ, J =3,4, k^O, 



(4.13) 



where, 



is a frequency cut-off operator. 

Knowing that exp (^— itAj {^^J—A^i, k^) are also L^— isometries, and using the fact that \E'(|^|) is 
L^-multiplier, we conclude, by interpolation. 



exp -itAj J-Ah, k ^ J -Ah [v 



< t 



- ~''c{^,P)\MLp'iR^), (4.14) 



for p > 2, i + — = 1, and for j = 1, 2, k e Z, j = 3,4, k 0, 
p p' 



4.3 Initial data decomposition 

In order to exploit the dispersive estimates derived in the preceding part, we have to find a suitable 
mollification of the initial data g^^^ and uq. To this end we consider a family of smooth functions 



il^5 e C^iO, oo),0<ips<l, ^5 / 1 as 5 ^ 0, 



and 



<p5 = Mxh) e C^{R^), < 05 < 1, 05 / 1 as 5 ^ 0. 
Finally, we regularize the data ^Jq^'', uq taking 

1 



^ |fc|<l/(5 



exp (— ifcxs) 



(4.15) 
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and, similarly, 

[uoj]s{xh,X3) = 4= 51 -^f-ixh [i's{\^\){uo,j(f)s){^,k)^ exp{-ikx3), j = 1,2,3. 

v2 \k\<l/5 



(4.16) 



In other words, we first multiply the data by a cut-off function to ensure integrability and then 
perform a similar cut-off in the frequency variable to ensure smoothness. We remark that 

• the functions ips are obviously multipliers for any 1 < p < oo, 

• the orthogonal projection V onto the kernel N'{B) commutes with the frequency cut-off repre- 
sented by ipSi 

• the operator iB represented in the Fourier variables by the matrix A commutes with the fre- 
quency smoothing, in particular, the time evolution of the mollified data remains restricted to 
the domain of frequencies bounded above as well as below away from zero. 



Finally, we write the initial data in the form 
= So,5 + go,5, where - ^hQo^s + go,<5 



curU [[uo]^]^ dx3 + 



[uo]^ = Vo,5 + vo,5, with [vq^s\i = -da;^qo,5, [vo,sh = dxiqo,s, 
and denote s^^s, ^e,5 the unique solution of the acoustic system 

£dtSs,s + diVa^-V^^^ = 0, 

supplemented with the initial data 
We have 

( 



(4.17) 
(4.18) 

(4.19) 
(4.20) 

(4.21) 



Q*(e,fc)exp 



t 

-i- 

e 



V 



\i{^,k) 

-Ai(^,A;) 

A3(e,A;) 

-A3(^,A;) 



where, in agreement with the previous observations, (i) each component of the vector 

Q(e,A;)te, V^,s\i^,k) 
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is of the form k), with ijj G C^(0, oo) for each fixed k e Z, (ii) as [so,5, Vo,^] G M[B]-^, 

^i(iei)«i(0 





Consequently, interpolating the decay estimates fl4.14l) with fl4.8p we may infer that 
Se,5 ^ in LP(0,T; W"''°°{Q)),Vs,5 in LP{0,T; W"''^{Q)) as e ^ 
for any fixed 6 > 0, any 1 < p < oo, and m = 0,1, . . .. 

5 Convergence, part II 



(4.22) 



We finish the proof of Theorem 12. II by means of another application of the relative entropy inequality 
fl3.2D. this time for the choice 



(5.1) 



where [se^5,V£,5] are solutions of the acoustic system (14.191 - I4.2ip . the properties of which were 
discussed in the previous section, while [g^, v^] is the solution of the target problem 



dt {Ahqs - qs) + V^qs ■ Vh {^hqs - qs) = 0, ^^(0, ■) = qo,s. 



(5.2) 
(5.3) 



cf. Proposition 12.11 



5.1 Initial data 



Going back to the relative entropy inequality (13. 2p with the "ansatz" (15. ip we get, in agreement with 
the hypotheses dm - EH]) , 

Se(go,e,no,e r(0, ■), U(0, ■)) (5.4) 



Jn2 



+ 



dx 



< C (||uo,e - [uo]5||i2(Q.R3) + ||^o!i - [Qi'']s\\L'i{Q.;Ri)^ 

c (||uo - [MsWhin.m + 11^0^^ - \Qi\\\h{n;Ri)) as £ ^ 0. 
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5.2 Viscosity 

We write 

r f §e{^xVe,5) : V^(U,,5 - U,) dx dt < Ci T f /i^ |V^(U,,5 - U,)| dx dt, 

Jo Jn Jo Jn 

where, by virtue of Korn's inequality, 

r / f^e |V^(U,,5 - U,)| dxdt< \ r [ (§e(V,U,) - §e(V^U,,5) : {\Je,& - ^e) dx dt + Ca/i,. 

JO Jn 2, Jo Jn 

5.3 Forcing term 

Furthermore, in accordance with the hypothesis fl2.17p . the convergence estabhshed in fl3.12p . and 
the decay estimates (14.221) . we get 

/ / Qe^xG ■ (Ue s-Ue) dx dt [ [ W xG ■ (v^ - u) dx dt = as 5 0, 
Jo Jn ' Jo Jn 

where we have used fl3.16p . (15. 2p . 

Summing up the previous estimates with (15.41) . we can rewrite the relative entropy inequality in 
the form: 



^,|u, - U,,5|2 + ^ {H{g,) - H'in,s)i0e - re,s) - Hin,s))) (r, •)] dx (5.5) 

< hi{e, S) + r f Qe {dt\Je,S + U, ■ V,U,,5) ■ (U,,5 - U,) dx dt 

Jo Jn 

+ - [ [ X Ue) ■ (lJe,S - Ue) dx dt 

e Jo Jn 

ire,s - Qs)dtH'{r^^s) + ^xH'ir^^s) ■ {re,sVe,5 - Qe^^ie)] dx dt 

~hlo -Pif^e,5)) dx dt. 

Here and hereafter, we use the symbol hi,i = 1,2,... to denote a function of e, 6 enjoying the 
following properties 



1 r 



Jo Jn 



hi{e, 5) hi{S) for e 0, with hi{S) for 5 0. 
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5.4 Estimating the remaining terms 

To begin, let us recall our convention that 

g = p'{g)=H"(:g) = l. 

Furthermore, in the following discussion, we make a systematic use of the dispersive decay estimates 
(5221). 

Step 1: 

We have 

{re,5 - Qe)dtH'{re,s) + VxH'{re,5) ■ {re,5^e,5 - Qe^e) - div^.Ue,5 - P{re,5)) (5.6) 

= pire,5) - P'{re.s) {re,5 - Qe) - div^Ue,5 

where, in accordance with fl4.19p . (15. 2p . 

dtre,5 + <^'^'^x{re,s^e,s) = ^dtQs + sdtS^^s + (iivx{re,s(ys + ^e,5)) = £dtqs + ediv^ {qs + Se,s)Ve,s ■ (5.7) 
Furthermore, we get 

QeV.H'ir) ■ (U,,5 - u,) = geV,[H\r,,s) - H"{g){re,s -q)- H'{g)] ■ (U,,5 - u,) (5.8) 

+ege^^V^qs ■ (U,,5 - u,) + ege^^V,Se,s ■ (U,,^ - u,) 
g g 

= geV.[H'{r,,s) - H"{l){r,^s - 1) - H'{1)] ■ {IJ.^s - u,) 

+ege^^Vxqs ■ {Ve,s - u^) - e'^ge(lJe,5 - Ue) ■ dtVe,5 
Q 



-ege{Ve,5 - Vie) ■ (w X V, 



Consequently, after a straightforward manipulation, the inequality (15.51) can be rewritten as 
follows: 

< hi{e, 6)+ [ [ ge (dt^s + ■ '^xVe,5) ■ (^6,8 - Ue) dx dt 
Jo Jn 



dx 
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/ / Qei^^ X V^) ■ (lJe,S - dx dt 

e Jo Jn 



+ - 



1 r 



Jo Jn 



p{rs,5) - p'irs,5){re,s - Qe) - P{Qe) 



divT-Up A dx dt 



1 



dx dt 

e Jo Jn 

0sV,[H'{r,^s) - H"{l){r,^s - 1) - H'{1)] ■ (U,,^ - u,) dx dt 

+- r I Qe^V.qs ■ (U,,5 - u,) dx dt 
e Jo Jn g 



Jo Jn 



Finally, we use the relation (15. 2 p to conclude that 
"1 , 1 



< hi{e, 6) + Qe {dt^s + Us ■ Vx-Ue,^) ■ {\Je,5 - dx dt 

Jo Jn 



dx 



1 r 



Jo Jn 

1 r 



P{re,5) - p'{rs^s){re,S - Qe) - P{Qe) 



divAJp I) dx dt 



e Jo Jn 



re,S - Qe)H"{re^s) dtQs + div^; ( (g^ + Se,5)V, 



dx dt 



Jo JQ 



QeVjH'{r,,s) - H"{l){r,,s - 1) - ^^'(1)1 ■ (U,,^ - u,) dx dt. 



Step 2: 

We rewrite the integral 



(5.9) 



f f Qe (dtvs + Ue ■ ^x^e.s) ' (^e.s - u^) dx dt = [ [ {df^s + ■ V^jV^) ■ (U^.^ - u^) dx dt 
Jo Jn ' ' Jo Jn 



+ / / Qe^S ■ ^x^e^S ■ (Ue,5 - U^) dx dt + / / Qe^ e,5 ' ^x^e^S ' {^e,S - U^) dx dt 

JO Jn Jo Jn 

Qe{^e,5 - ^e) ' V^rUe^^ ■ {\J e,5 - ^e) dx dt, 



Jn 



where, by virtue of the uniform bounds (13.11 - 13^ . combined with the dispersive estimates (I4.22p . 

/ / Qe^S ' V X ^e,5 ■ (U£,5 - Ug) dx dt + / / QeVe,& " '^x^e,& " {^e,6 - U^) dx dt| -> aS £ ^ 

10 Jn Jo Jn 
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uniformly in r G [0,T]. 

Consequently, we may infer that 



+ / Qe {dt^s + V5 ■ V^v^) ■ {JJs^s - Uj) dx dt 
Jo Jn 



+- 



1 r 



Jo Jn 



P{re,5) - P'{r^^s)ire,5 - Qe) - P{.Qe) 



div,,U=. A dx dt 



£ JO Jn 

1 r 



dx dt 



Jo Jn 



QsV, H'{r,^s) - H"{l){r,,s - 1) - i^'(l) ■ {Ue,s - u,) dx dt. 



Step 3: 

In view of the uniform bounds ( 13.41 - I3.6P , we have 



dx 

(5.10) 



ess sup — \\[p{rs,5) - p'{r,^s){re,s - &) - p{Qe)] (tr^LUn) ^ 
te(o,T) £ ^ ' 



whence, by virtue of the dispersive estimates fl4.22p . 



Jo Jn 



p{re,5) - p'{re,s){re,s - 0s) -p{Qe) div^.U£,5 dx dt 



Jo Jn 

uniformly for r G [0, T]. 
Similarly, we obtain 



P{.re,5) - P'{r^^s){re,5 - Qe) -p{.Qe) 



div^.Ve A dx dt — )■ as £ — )■ 0, 



H\n,s)~H\l){n^,-l)-H'{l) 



(5.11) 



= - {H"{r,^s) - H"{1)) V,{qs + s,,s) ^ H'"{l)qsV,qs in L°°(0, T; n L°^(l]; i?=^)) as £ ^ 0. 
Thus, using f l3.12p . f l3.13p we conclude that 



H'ir.^s) - H"{l){r,^s - 1) - ^^'(1) ■ (U,,^ - u,) dx dt -> as £ ^ (5.12) 



JO Jn 
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uniformly in r G [0,T]. 

In view of flS.lip . f l5.12p . the relative entropy inequality (l5.1Up reduces to 



< 



/l3(£:,(5) + SUPig[oy]||v5(t, ■)||h/i,oo(^2.^2) f Se(^ge,Uere,5,^e,s) {t, ■) dt 

•J 



da; 

(5.13) 



+ / j Qe (dtVS + V5 ■ Wx^^s) ■ (^6,5 - Ue) dx dt 

Jo Jn 



1 r 



e Jo Jn 



dx dt. 



step 4: 

Similarly to the above, we deduce from fl5.13p that 

■1 , ,0 1 



J ^ _ ^ s 



dx 

(5.14) 



< h4{e,5) + SUPtg[o,r]||v5(t, ■)||^yi,oo(^2.^2) Se [Qs^-^s ^^£,5, U^^^j (t, ■) dt 

{dt^s + ■ Vx^s) ■ (v5 - u) dx dt + ^ ^(g^ - g^^^) dtqs + dwxiqs^s) 



dx dt, 



where 



{^t^rs + V5 ■ V^v^) ■ (v^ - u) dx + l {qs - g'^^'^) dtqs + div^ 
1 d 



qs^s 



dx 



- ^ (v5 • VxV^ ■ u + g^'^^ dXN x{q8^ s)) dx. 

Furthermore, we have 

div^(g5V5) = \/xqi ■ V5 = -Vxq& ■ V^g5 = 0, 
and, by virtue of fl2llD . (ISlij) . 

- / (9tV5 ■ u + Stg^^'^^O dx = dt {Ahqs - qs) g^^^ dx = - / g'^^^VhiAhqs) ■ dx 
Jn ^ ^ Jn Jn 

= - {u X u) ■ vs^hqs dx. 
Jn 
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On the other hand, a routine manipulation gives rise to 

V5 ■ Vx-V5 ■ u + X u) ■ \sAhqs = u ■ V^lv^l^; 

whence 

{dt^s + V5 ■ V^v^) ■ (v5 - u) dx+ / {qs - g'^^'^) \dtqs + div^fg^v^ 
n Jn 



dx 



Finally, multiplying the equation (15 ■3p on we obtain 

t;^ f (|V/ig5p + qs) dx = [ V ■ Wh{Ahqs)qs dx = - [ v ■ Vhqs^hqs dx = 0. 
2 dr ^ ^ Jn Jn 

In view of the previous discussion, the relation (15 .14^ reduces to 

/" r 1 1 

Se(^ge,Usr^^s,Us,s) = -ge\ue - IJe^'^ + ^[H^g^) - H'{re,s){ge - re,5) - H{ 



< h^{e,6) + snpt^[Q^T^\\vs{t,-)\\whoo{R^;R^) / SJge,u,r,^s,'Ue,5) (t,-) dt 



with 



dx 

(5.15) 



hii^e, 6) /i4(5) for s 0, with /i4(5) for 5 0. 

Passing to the limit, first for e — )■ and then for 5 — )■ 0, in (15.151) . we complete the proof of 
Theorem 12.11 
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